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摘 要 : 利用 新 疆 地 区 2016 年 夏季 (6 一 8 月 ) 近 2000 个 


276400) 


区 域 自动 站 小 时 降水 数据 ,以 最 新 一 代 GPM(Global Precipi- 


tation Measurement Mission , 全 球 降水 测量 计划 )IMERG(Integrated Multi-satellite Retrievals for GCPM) 卫 星 降 水 产品 为 


初始 场 , 运 


概率 密度 匹配 与 最 优 插值 两 步 融 合 校 正方 法 (probability density function-optimal interpolation, PDF- 


01) ,开展 干旱 区 多 源 降水 融合 试验 。 首 先 ,通过 将 实况 观测 与 卫星 反 演 降水 进行 概率 密度 函数 (PDF) 匹 配 ,从 而 实 
现 对 IMERG 降水 的 偏差 订正 ;然后 以 偏差 订正 后 的 IMERG 降水 为 初 佑 场 , 实 况 降水 为 观测 场 ,通过 最 优 插 值 (01) 


将 偏差 订正 后 的 IMERG 降水 量 估 计 值 与 雨量 计 实 况 相 


二 
结合 


,最 终 得 到 经 过 两 步 校 正 后 的 新 疆 地 区 小 时 雨量 数据 。 


交叉 验证 结果 表明 ,与 雨量 计 实 况 和 原始 卫星 降水 数据 相 比 ,PDF-OI 两 步 校正 方法 构建 的 融合 降水 不 仅 大 大 消除 


了 系统 误差 ,而 且 显 车 提高 了 数据 精度 。 
关键 词 : 星 - 地 联合 ; 小 时 降水 ; 融合 ; 新 疆 


降水 是 与 大 气 潜 热 释放 有 关 的 大 气 活动 , 降 水 
过 程 使 得 地 球 表 面 接收 的 太阳 辐射 被 重新 分 配 到 
大 气 层 内 部 ,从 而 驱动 大 尺度 环流 "。 掌 握 一 个 特 
定 区 域 的 降水 时 空 分 布 有 助 于 更 好 地 理解 这 一 区 
域 的 生态 .环境 变化 ,而 这 其 中 最 关键 的 是 准确 获 
取 整 个 区 域 的 降水 数据 。 由 于 降水 在 不 同 空间 和 
时 间 上 有 很 大 的 差异 ,获得 精确 的 区 域 和 全 球 降水 
言 息 ,往往 是 一 项 艰巨 的 挑战 。 

当前 ,获取 降水 信息 的 手段 主要 有 3 种 : 台 站 实 
况 观测 ”雷达 估 测 和 卫星 遥感 反 演 。 站 点 观测 具 
有 较 高 的 精度 ,但 台 站 大 多 位 于 低 海 拔 地 带 , 且 分 
布 不 均 ,不 能 有 效 反 映 降 水 的 空间 变化 特征 ,空间 
代表 性 不 强 ;基于 气象 雷达 回 波 推算 的 降水 数据 具 
有 一 定 区 域 佑 测 的 优点 ,然而 Z-[( 反 射 率 - 雨 强 ) 关 
系 在 不 同 的 降水 系统 .季节 、 区 域 都 有 很 大 差异 , 影 
啊 了 雷达 佑 测 降 水 的 准确 性 ,此 外 ,雷达 探测 易 受 
到 地 物 阻挡 波束 效应 以 及 速度 模糊 等 限制 而 使 降 
水 佑 测 精度 不 高 ,实际 应 用 中 仍然 存在 较 大 限制 ; 
遥感 和 地 理 信 息 科 学 的 快速 发 展 , 为 大 范围 降水 同 
步 观测 提供 了 新 方法 ,卫星 反 演 降水 产品 具有 全 天 
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候 全 球 履 盖 以 及 准确 反映 降水 空间 分 布 的 独特 优 
BASS! ,使 其 成 为 系统 了 解 区 域 乃 至 全 球 降水 情况 及 
其 变化 的 重要 手段 。 

卫星 遥感 具有 大 范围 空间 连续 探测 的 优势 ,但 
由 于 其 是 非 直接 观测 ,卫星 反 演 降水 的 物理 原理 和 
算法 还 存在 很 多 局 限 性 ,使 得 卫星 降水 产品 具有 很 
大 不 确定 性 ,精度 相对 较 低 。 因 此 ,在 应 用 之 前 ,有 
必要 进行 卫星 降水 产品 的 订正 。 针 对 卫星 降水 产品 
精度 较 低 ,许多 研究 致力 于 发 展 卫星 降水 与 地 面 观 
测 实况 的 融合 方法 ,其 中 包括 客观 统计 分 析 “ ,最 优 
插值 技术 ”中 ,神经 网 络 呈 号 ,双核 平滑 "和 地 理 差异 
分 析 “ 等 。 其 中 ,国家 气象 信息 中 心 针 对 多 源 降水 
融合 开展 了 大 量 人 研究 工作 "1, 由 其 研制 的 降水 融合 
产品 (CMPA) 已 通过 全 国 综合 气象 信息 共享 平台 下 
发 至 各 级 气象 部 门 ,在 业务 及 科研 中 均 得 到 广泛 应 
FA), RT , CMPA 融合 降水 主要 以 Climate Predic- 
tion Center (CPC) Morphing Technique (CMORPH)”” E 
星 降 水 产品 为 初始 场 , 通 过 对 该 产品 在 新 疆 地 区 与 
其 他 卫星 降水 产品 进行 对 比 评估 ,结果 显示 最 新 一 
代 全 球 卫 星 降 水 产品 IMERG 在 干旱 区 的 新 疆 比 
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CMORPH 精度 更 高 ,尤其 在 小 时 尺度 ”。 因 此 ,本 研 
究 在 CMPA 融合 降水 的 研究 方法 基础 上 ,在 新 疆 地 
区 以 IMERG 降水 为 初始 场 开 展 多 源 降水 融合 试验 ， 
以 期 获得 更 高 质量 的 新 疆 地 区 小 时 降水 融合 产品 。 

开展 新 疆 地 区 卫星 降水 产品 与 地 面 观 测 实况 
的 融合 试验 ,对 准确 反映 新 疆 降水 分 布 信息 和 降水 
演变 过 程 特征 意义 深远 。 通 过 研制 的 高 时 空 分 辨 
率 格 点 降水 数据 集 , 将 有 助 于 推进 我 们 对 新 疆 地 区 
小 时 强 降 水 特征 和 发 生发 展 规律 的 了 解 ,为 该 地 区 
提高 预报 准确 率 、 预 防 气象 次 生 灾害 提供 数据 支 
撑 , 同 时 ,为 今后 研制 中 亚 区 域 高 时 空 分 辩 率 降水 
产品 奠定 技术 基础 。 


1 数据 的 获取 


1.1 GPM IMERG 卫星 降水 
GPM(Global Precipitation Measurement Mission , 


计 比 TRMM 卫星 多 了 4 个 高 频 波段 ,从 而 对 小 十 和 
降雪 更 加 敏感 。 与 主要 观测 中 、 大 型 降水 的 TRMM 
卫星 相 比 ,GPM 提 高 了 小 雨 和 固态 降水 的 观测 能 
而 这 2 种 是 新 疆 地 区 的 主要 降水 类 型 ,因此 ,对 本 人 研 
究 区 来 说 ,GPM 的 观测 成 果 有 具有 重要 人 研究 意义 。 

本 人 研究 选取 GPM 的 3 级 IMERG 产品 为 降水 初 
始 场 ,时 空 分 辨 率 为 0.5 bh、0.1°x0.1°。 研 究 中 将 0.5 
h 降 水 累加 为 1h 与 逐 小 时 实况 降水 匹配 ,时 间 段 为 
2016 年 6 一 8 月 。 

1.2 实况 降水 

新 疆 地 区 近 2000 个 区 域 自 动 站 逐 时 降水 数据 ， 
实况 降水 (命名 为 OBS) 时 间 范 围 与 IMERG 降水 数 
据 一 致 。 实 况 数 据 由 新 疆 气 象 局 信息 中 心 整理 并 
通过 气候 极 值 检验 、 单 站 极 值 检验 和 数据 一 致 性 检 
验 等 质量 控制 。 为 保证 建 模 和 验证 的 独立 有 效 性 ， 
E AIR ir SE He DX. 19 个 国际 交换 站 。 由 于 新 疆 区 域 
自动 站 大 部 分 为 翻斗 式 雨量 计 , 只 能 测量 降雨 ,在 


全 球 降 水 测量 计划 ) 是 由 美国 国家 航空 航天 局 
(NASA) 与 日 本 宇宙 航空 研究 开发 机 构 (JAXA ) 联 合 
开发 的 全 球 降水 测量 系统 ,是 TRMM(Tropical Rain- 
fall Measuring Mission ,热带 降雨 测量 任务 ) 的 后 续 卫 
星 降水 计划 。GPM 是 由 1 颗 主 卫星 和 10 颗 小 卫星 
组 成 的 1 个 卫星 群 ,每 3h 获 得 1 次 全 球 降水 分 布 数 
据 资 料 。GPM 核心 卫星 搭载 了 双 频 雷达 和 微波 辐 
射 计 ,其 中 , 双 频 雷达 比 TRMM 卫星 多 出 一 个 Ka 波 
段 ,可 探测 更 为 详细 的 降水 微 物理 信息 ;微波 辐射 
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冷 季 停止 观测 , 故 研 究 选 取 暧 季 的 6 一 8 月 开展 研 
究 。 图 1 显示 了 新 疆 区 域 0.1°x0.1° 网 格 内 台 站 数量 
分 布 情况 ,新 疆 地 区 0.1° 网 格 共有 18038 个 ,其 中 
1509 个 为 有 效 格 点 (<10%), 即 该 格 点 内 至 少 包 含 1 
个 雨量 站 。 对 1509 个 有 效 格 点 进一步 细 分 ,包含 1 
个 雨量 站 的 有 1267 个 ( 占 84%),2 个 雨量 站 及 2 个 
以 上 的 分 别 有 199 个 和 43 个 , 绝 大 部 分 有 效 格 点 都 
只 含 1 个 雨量 站 。 按 照 10 折 交叉 验证 ,将 有 效 格 点 
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图 1 新 疆 区 域 0.1? x 0.1° 网 格 台 站 个 数 分 布 示 意图 


Fig. 1 Distribution of gauge numbers in the 0.1° x 0.1° lat/lon grid boxes in Xinjiang 
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随机 分 为 10 组 ,每 次 选取 9 组 用 于 PDF 订正 和 0I 融 
合 , 剩 余 的 1 组 作为 独立 数据 集 进行 偏差 订正 后 以 
及 最 终 融合 产品 的 精度 验证 ,如 此 重复 10 次 ,使 每 
个 台 站 都 独立 地 进行 了 校正 和 检验 ,以 此 保证 训练 
样本 及 验证 样本 的 代表 性 和 客观 性 。 


2 研究 方法 


2.1 建立 有 效 格 点 数据 对 

对 新 疆 区 域 1509 个 有 效 格 点 进行 卫星 数据 与 
地 面 实况 的 时 空 匹 配 ,将 落 入 同一 网 格 的 所 有 实况 
站 点 降水 取 平 均 与 该 网 格 卫 星 降水 形成 1 个 空间 匹 
配 数 据 对 ,时 间 上 则 取 对 应 的 同一 时 次 进行 匹配 ， 
处 理 过 程 中 需 注意 卫星 降水 为 世界 时 , 需 加 8 与 
实况 降水 的 北京 时 对 应 ,将 配对 好 的 数据 样本 作为 
下 一 步 PDF(probability density function) 订 正 与 0I 
(optimal interpolatiom) 融 合 的 数据 基础 。 

2.2 PDF 匹配 法 误差 订正 

PDF 订 正 主要 是 以 雨量 计 观 测 的 实况 降水 概率 
密度 对 卫星 降水 佑 计 的 概率 密度 进行 “标定 ”, 使 两 
者 的 概率 密度 值 相 匹配 ,从 而 达到 去 除 卫 星 降水 系 
统 误差 的 目的 。 具 体 步 又 如 下 : 

(1) 由 于 小 时 降水 具有 零 降水 发 生 概 率 大 以 及 
更 小 时 空 尺度 等 特点 ,为 保证 获得 稳定 的 累积 概率 
密度 ,以 待 订正 的 IMERG 降水 格 点 (以 下 称 当前 格 
点 ) 为 中 心 ,4?x4? 空 间 范 围 为 空间 窗口 选取 有 效 格 
点 数据 对 ,同时 ,针对 逐 时 降水 的 日 变化 特征 ,以 当 
前 时 间 为 基准 , 取 当 前 日 期 向 前 20 d 以 及 当前 时 刻 
向 前 6 个 时 次 为 时 间 和 窗口 , 即 每 个 订正 时 次 共 取 当 
前 格 点 4?x4? 空 间 范 围 内 有 效 格 点 的 120 个 匹配 时 
次 作为 订正 样本 ,从 而 收集 得 到 匹配 的 地 面 观测 和 
卫星 降水 进行 PDF 订正 ,详细 过 程 可 参考 文献 [9]。 

(2) 在 上 述 得 到 有 效 格 点 数据 对 基础 上 ,选取 
实况 观测 和 卫星 降水 均 为 非 缺 测 的 作为 有 效 数 据 
样本 ,针对 每 个 样本 ,以 到 当前 时 刻 .当前 格 点 的 时 
空 距离 为 权重 进行 最 有 效 样本 的 筛选 。 

(3) 按 同一 降水 密度 等 级 ,分 别 计算 每 个 当前 
格 点 时 空 匹配 的 有 效 数据 样本 中 实况 降水 与 卫星 
降水 的 累积 概率 密度 ,由 此 得 到 实况 降水 与 卫星 降 
水 的 累积 概率 密度 匹配 表 , 进 一 步 计 算 当 前 格 点 卫 
星 降水 值 在 卫星 降水 累积 概率 密度 表 中 所 处 的 密 
度 等 级 ,得 到 相应 的 卫星 降水 概率 密度 值 ,同时 , 计 
算 该 卫星 降水 概率 密度 值 所 匹配 的 地 面 降水 累积 


概率 密度 及 对 应 的 实况 降水 量 ( R,), 则 订正 后 的 卫 
星 降 水 值 即 为 R,。 
2.3 最 优 插值 (OI) 

最 优 插 值 分 析 以 偏差 订正 后 的 IMERG 降水 为 
初 佑 场 ,以 站 点 实况 降水 为 真 值 ,每 个 格 点 上 的 最 
终 降 水 分 析 值 4 由 两 部 分 组 成 ,一 部 分 为 该 点 的 初 
估 值 有 尺 , 即 卫星 反 演 降水 量 , 另 一 部 分 为 该 格 点 卫 
星 降水 值 与 实况 降水 量 的 偏差 ,由 于 大 部 分 网 格 内 
都 没有 实况 站 点 ,因此 ,这 个 偏差 由 该 格 点 周 于 一 
定 范围 内 nn 个 有 效 格 点 数据 对 的 偏差 加 权 估 计 得 
到 ,公式 如 下 : 


A,=F,+¥W(0,-F)) (1) 
i=1 


式 中 法 和 ;分 别 为 当前 分 析 格 点 和 周围 有 效 格 点 ; 
唉 为 有 效 格 点 数据 对 的 权重 函数 ,表示 ;点 上 实况 
值 与 初 估 值 的 偏差 在 订正 时 分 配 的 权重 。 值 得 注 
意 的 是 由 于 新 疆 站 网 稀 琉 ,将 分 析 半 径 设置 为 100 
km ,以 保证 在 本 区 域 能 搜索 到 一 定数 量 的 有 效 格 
点 ,再 以 距离 由 近 及 远 排 序 , 从 中 选取 距离 分 析 格 
点 最 近 的 个 有 效 格 点 参与 最 优 插值 ,n 取 值 为 9。 
权重 系数 取 由 该 分 析 点 降水 值 4 的 最 小 误差 方差 
决定 。 

假设 观测 场 与 初 估 场 均 是 无 偏 且 观 测 误差 和 
初 估 场 误差 不 相关 "", 则 式 (1) 的 权重 系数 取 可 以 
通过 求解 下 面 的 线性 方程 组 来 得 到 
fr- (2) 


J 


SO i, [ENA S , AL isj=1, 2……… ,n; uf 
几 分 别 表示 初 估 场 误差 协 相 关 和 实况 观测 误差 协 
相关 ; A; CA) ) A, 为 ;点 0 点) 观测 误差 标准 差 wy ( 0?) 
和 初 估 场 误差 标准 差 o! (o ) 的 比率 。 在 0I 融 合算 
法 中 ,关键 是 量化 初 估 场 和 观测 场 的 误差 结构 ,这 
就 要 求 式 (2) 中 的 必 ae. o Alo! 均 为 已 知 量 。 
本 研究 通过 对 2016 年 6 一 8 月 的 样本 数据 进行 统计 
得 到 上 述 参 数值 。 参 数 统计 过 程 中 计算 得 到 的 卫 
星 反 演 降 水 (IMERG ) 的 均 方 误差 与 雨 强 的 散 点 分 
布 如 图 2 所 示 , 将 有 效 格 点 中 任意 2 格 点 卫星 估 测 
值 与 实况 值 的 误差 与 距离 做 相关 分 析 , 可 得 到 初 估 
场 的 误差 相关 曲线 (图 3)。 

上 述 参 数 中 , yi 由 图 3 拟 合 
定 初 估 场 2 格 点 i 和 j 的 距离 则 通 


函数 式 得 到 , 即 给 
过 图 3 可 得 到 2 格 
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图 2 不 同 雨 强 下 卫星 反 演 降水 (IMERG ) 的 均 方 误差 
Fig. 2 Mean square errors of the IMERG under different 


rainfall intensities 


oo [ y=A1 xexp(—x/t1)+y0 
pel A1=0.699 
: y0=0.0086 
11=28.55 
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图 3 任意 2 格 点 误差 相关 
Fig. 3 Error correlation for the bias-corrected IMERG as a 


function of separation between two grid boxes 


点 间 的 误差 协 相关 ; uo 为 实况 观测 误差 协 相 关 ,前 
提 是 每 个 台 站 的 观测 误差 只 与 自身 的 观测 有 关 , 由 
此 可 得 出 ,iy 相 同时 jw 为 1, 不 同时 pi 为 0; of R 
示 初 估 场 误差 标准 差 , 可 由 图 2 拟 合 函 数 通过 给 定 
卫星 降水 值得 到 , 同 理 o. 可 通过 实况 观测 的 均 方 误 
差 与 南 强 的 散 点 分 布 计算 得 到 。 详 细 计 算 过 程 可 
参考 文献 [7]。 

通过 求解 式 (2) 确 定 权 重 系数 WW, 最 后 由 式 (1) 
得 到 最 终 降 水 分 析 值 41。 
2.4 检验 评估 

台 站 实测 数据 是 最 真实 的 降水 监测 ,本 研究 按 
10 折 交叉 验证 利用 独立 站 点 进行 融合 降水 产品 的 
效果 检验 ,并 与 原始 IMERG 小 时 降水 的 评估 结果 进 
行 对 比分 析 。 主 要 选取 4 个 统计 指标 进行 评估 ,分 
为 两 类 ,第 一 类 相关 系数 (Correlation Coefficient, 
CC) , 它 反映 融合 降水 与 实测 降水 的 一 致 性 程度 ;第 
二 类 绝对 平均 误差 (Mean Error, ME) 、 相 对 误差 


(Relative Bias, RB) 、 均 方 根 误差 (Root Mean 
Squared Error, RMSE) ,用 来 描述 融合 降水 对 比 实况 
的 误差 大 小 。 公 式 如 下 所 示 : 


n 


DL: -#)(y;-7) 


cC=— - (3) 
Se- $o- 
ME = 1E («-7) (4) 
; (x, -y;) 
RB= =" x 100% (5) 
2 
RMSE = 1Y>(% -y (6) 


式 中 :CC 为 相关 系数 ;ME 为 绝对 误差 ;RB 为 相对 误 
差 ;RMSE 为 均 方 根 误差 ，#=》 ，7= TD ys n 


为 样本 容量 ; x, y, 分 别 为 降水 估 测 值 和 气象 站 观 
测 值 。 


3 ARAM 


3.1 PDF 匹配 对 原始 IMERG 降水 的 改进 

通过 与 实况 降水 概率 密度 匹配 ,对 原始 IMERG 
卫星 降水 进行 偏差 订正 。 图 4 显示 交叉 验证 后 的 实 
况 (OBS)、 原 始 IMERG、 以 及 偏差 订正 后 的 降水 ( 命 
名 为 pdf_CRT)3 种 降水 数据 的 概率 密度 分 布 ,注意 
由 于 小 时 降水 中 0 值 降水 占 绝 大 多 数 比 例 , 因 此 ,图 
4 从 大 于 0 mm 降水 开始 统计 。 图 4 中 反映 出 原始 


20 maOBS CIMERG (Copdf CRT 


‘dade 


(0~0.1] 0. 1~0.2](0.2~0.4(0.4~0.6] (0.6~1] (1~2] (2~5] >5 
降水 量 /mm 


概率 密度 分 布 /% 


图 4 2016 年 夏季 新 疆 地 区 实况 (OBS) .订正 前 (IMERG)、 
订正 后 (pdf CRT) 小 时 降水 概率 密度 分 布 对 比 
Fig. 4 Comparison of the PDF's of hourly precipitation among 
OBS, original IMERG, and bias-corrected 
IMERG in summer 2016 over Xinjiang 
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IMERG 降水 在 大 多 数 情况 下 的 概率 分 布 均 远 远 高 
于 实况 ,只 在 降水 量 较 大 时 (>2 mm.h) 低 于 实况 降 
水 ,而 经 过 概率 密度 匹配 订正 后 的 pdf_CRT, 在 各 个 
降水 等 级 的 概率 密度 值 都 明显 与 实况 更 为 接近 , 具 
体 表现 为 在 原始 IMERG 的 PDF 偏 大 时 减 小 其 概率 
密度 值 ,而 在 原始 IMERG 的 PDF 偏 小 时 增 大 其 概率 
密度 值 ,体现 了 PDF 订 正方 法 中 实况 降水 概率 分 布 
对 卫星 降水 的 调节 作用 。 表 1 的 定量 评估 结果 也 显 
ZR , PDF 订正 后 相对 误差 由 订正 前 的 48.41% 减 小 到 
23.77% ,降幅 达 50%。 表 1 还 显示 出 PDF 订正 后 相 
关系 数 比 订正 前 略 有 下 降 ( 下 降 0.006) ,表明 PDF 订 


1.0 


0.8 


0.6 


相关 系数 


0.4 


0.2 


正方 法 主要 起 到 了 降低 卫星 降水 系统 偏差 的 作用 ; 
进一步 给 出 订正 后 降水 数据 的 偏差 随时 间 变 化 的 
曲线 (图 5b ~5d), 图 中 显示 pdf_CRT 降 水 的 ME 与 
RB 在 各 个 时 有 段 都 比 原始 IMERG 降水 显著 降低 , 系 


表 1 2016 年 6 一 8 月 3 种 降水 产品 交叉 验证 结果 


Tab. 1 Results of cross-validation of three precipitation 


products from June to August 2016 


相关 ”平均 误差 相对 误差 ” 均 方 根 误差 
系数 /(mm"h"') /% /(mm-h"') 
IMERG 0.288 0.15 48.41 1.20 
pdf_CRT 0.282 0.08 23.77 1.24 
pdf-oi_CRT 0.379 0.07 21.43 1.18 
pdf CRT ——pdfoi CRT 


平均 误差 /mm 
= 


-0.3 | 14 
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06-10 06-20 06-30 07-10 07-20 07-30 08-09 08-19 08-29 
pdf CRT 


—— pdf-oi CRT 


1 


©  ....- IMERG 


相对 误差 /% 
S 
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图 5 2016 年 6 一 8 月 独立 检验 统计 值 的 时 间 变 化 


Fig. 3 Time series of (a) CC, (b) ME, (c) RB and (d) RMSE between the gauge-based analysis and the IMERG, pdf_CRT, 


and pdf-oi_CRT precipitation over Xinjiang from June to August 2016 
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统 偏差 更 接近 0。 
3.2 OI 的 改进 效果 

由 表 1 可知 ,在 各 项 指标 中 ,经 过 2 步 校正 的 降 
水 产品 (命名 为 pdf-oi_CRT) 都 表现 最 优 , 与 原始 
IMERG 相 比 ,相关 系数 提高 31.6% ,绝对 平均 误差 下 
降 53.3% ,相对 误差 下 降 55.7% , 均 反 映 出 大 幅 改 进 
效果 。 由 图 $ 可 知 ,3 种 降水 产品 各 种 评估 指标 的 时 
间 变 化 曲线 ,直观 地 反映 出 pdf-oi_CRT 降 水 质量 在 
绝 大 部 分 时 间 点 均 比 原始 IMERG 降水 有 很 大 提高 。 
如 CC 从 订正 前 的 0.06~ 0.56 提 高 到 PDF-OLI 订 正 后 
的 0.06 ~ 0.66; ME FI RB 分 别 由 订正 前 的 -0.19 ~ 
0.42 mm +h’! #il-27.14% ~ 658.22% 下 降 到 -0.14~ 
0.27 mmh” 和 -23.06% ~ 413.2%. 

以 上 分 析 均 是 基于 单个 统计 指标 的 对 比 评 


标准 差 /mm 


准 差 /mm 


标 


1.0 
标准 差 /mm 


佑 ,图 6 给 出 了 各 种 降水 产品 与 实况 的 泰勒 图 , 泰 
勒 图 2 以 网 形 化 的 方式 提供 了 1 种 针对 模式 (或 估 
测 ) 精 度 进 行 综合 评价 的 方法 。 通 过 模式 产品 之 间 
的 相关 系数 . 均 方 根 误差 和 标准 差 来 量化 它们 之 间 
的 相似 性 , 离 实况 (字母 A) 距 离 越 近 则 表明 一 致 性 
越 好 。 泰 勒 图 对 评估 多 个 复杂 模型 或 测量 许多 不 
同 模型 的 整体 性 能 方面 特别 有 用 。 从 图 6 可 以 看 
出 ,结果 分 布 较 一 致 , 离 实况 最 近 的 均 是 pdf-oi_CRT 
(字母 D) ,表明 pdf-oi_CRT 降 水 的 评 佑 结果 最 好 。 
进一步 分 析 如 图 6a 所 示 ,在 整个 研究 时 段 内 实况 降 
水 的 标准 差 ( 黑 色 虚 线 ) 为 1.16 mm, 而 IMERG、 
pdf_CRT 和 pdf- oi CRT 分 别 为 0.74、0.86 和 0.92 
mm, 这 意味 着 pdf-oi_CRT 降水 的 标准 差 最 接近 实 
况 。 同 时 ,pdf-oi CRT 与 实况 降水 的 CC 最 高 


标准 差 /mm 
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5 


示 


hi 


0.5 1.0 1.5 2.0 
标准 差 /mm 


TE: A 表示 实况 降水 ,B 表 示 IMERG,C 表 示 pdf_CRT,D 表 示 pdf-oi_CRT。 
图 6 3 种 降水 产品 2016 年 整个 研究 时 段 以 及 6 一 8 月 的 泰勒 图 


Fig.6 Taylor diagrams of the three precipitation products for the whole period and each month from June to 


August 2016 over Xinjiang 
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(0.39) ,而 RMSE( 绿 色 虚 线 ) 均 低 于 其 他 2 种 产品 。 
同样 的 结果 也 可 以 从 各 月 的 泰勒 图 得 到 反映 (图 
6b ~ 6d). 
3.3 降水 个 例 分 析 

2016 年 7 月 31 日 至 8 月 1 日 新 疆 地 区 自 西向 东 
出 现 了 大 范围 降水 ,其 中 伊犁 地 区 “7.317 暴 雨 最 为 
显著 ,此 次 大 降水 天 气 过 程 由 7 月 31 日 17:00 开 始 
至 8 月 1 日 17:00 结 束 , 雨 强 在 8 月 1 日 02:00 左 右 达 
到 最 强 。 图 7a~7d 分 别 给 出 了 8 月 1 日 03:00( 北 京 
时 ) 的 地 面 自动 气象 站 观测 、IMERG 卫星 降水 、 
pdf_CRT 订 正 降 水 以 及 pdf-oi_CRT 融 合 降 水 分 布 。 
从 图 7 可 看 出 ,实况 降水 此 时 主要 集中 于 伊犁 E 
州 ` 塔 城 地 区 ,而 大 值 中 心 在 伊犁 地 区 ;图 7b IMERG 
卫星 反 滨 降水 显示 降水 区 域 远 远 大 于 实况 分 布 , 且 
降水 大 值 中 心 比 实况 偏 北 位 于 塔 城 地 区 托 里 县 一 
带 , 此 外 ,阿勒泰 .阿克苏 以 及 喀什 地 区 均 有 大 面积 
降水 ,而 实况 站 点 显示 该 区 域 仅 有 少量 台 站 出 现 降 
水 ,这 反映 出 卫星 反 演 降 水 普遍 存在 的 高 值 低估 、 
低 值 高 估 的 现象 ;图 7c 是 经 过 概率 密度 匹配 订正 后 
得 到 的 格 点 降水 pdf-CRT, 可 以 看 出 ,由 PDF 偏差 订 
正 后 ,原始 卫星 降水 的 虚报 和 高 估 现 象 得 到 显著 改 


TS°E 80°E 85°E 90°E 95°E 


善 ,上 且 降水 大 值 中 心 也 下 移 至 与 实况 一 致 的 伊犁 地 
区 ,然而 不 足 之 处 是 pdf-CRT 降 水 过 滤 掉 了 许多 实 
况 站 发 生 降水 的 区 域 ,尤其 在 伊犁 , 博 州 西部 地 区 ; 
图 7d 为 概率 密度 匹配 与 最 优 插值 两 步 校正 最 终 得 
到 的 pdf-oi_CRT 融 合 降水 ,该 图 显示 原始 IMERG 降 
水 在 PDF 偏差 校正 基础 上 再 融合 实况 降水 ,不仅 修 
正 了 图 7c 中 PDF 处理 后 被 错误 过 滤 的 有 降水 区 域 ， 
而 且 确 保 了 在 有 实况 站 点 的 区 域 与 实况 保持 一 致 ， 
在 无 实况 观测 区 域 保 留 卫星 反 演 降水 分 布 ,有 效 结 
合 了 2 种 观测 方式 的 优点 。 
3.4 同类 融合 降水 产品 对 比 

本 研究 是 建立 在 文献 [8-12] 的 方法 基础 上 ,不 
同 之 处 在 于 ,本 研究 选取 在 新 疆 地 区 精度 更 高 的 
GPM IMERG 卫星 降水 产品 作为 初始 场 ,日 在 算法 处 
理 过 程 中 根据 新 疆 地 区 降水 特点 进行 了 本 地 化 处 
理 , 因 此 ,有 必要 将 得 到 的 pdf-oi_CRT 融 合 降 水 与 
家 气象 信息 中 心 研发 的 全 国 范 围 降水 融合 产品 CM- 
PA 进行 比较 。 需 要 指出 的 是 ,由 于 缺乏 独立 验证 数 
据 集 ,无 法 将 pdf-oi_CRT 产 品 与 CMPA 产品 进行 全 
面 的 对 比 评估 ,本 研究 通过 将 两 者 在 空间 分 布 上 进 
行 对 比 , 分 析 评 估 两 者 的 差异 。 同 样 以 “7.31” 暴 雨 
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图 7 2016 年 8 月 1 日 03:00 新 疆 地 区 雨 强 分 布 


Fig.7 Hourly precipitation at 03:00 (UTC+8) 1 August 2016 over Xinjiang 
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8 2016 年 8 月 1 日 03:00 新 疆 地 区 雨 强 分 布 对 比 
Fig. 8 Comparison of the hourly precipitation at 03:00 (UTC+8) 1 August 2016 over Xinjiang 


过 程 为 例 , 图 8 给 出 融合 降水 CMPA( 图 8e) 及 其 初 
始 场 CMORPH( 图 8d) 同 期 降水 空间 分 布 ,与 本 文 初 
始 场 IMERG (图 8b) 以 及 融合 降水 pdf-oi_CRT( 图 
8c) 进 行 对 比 。 图 8d 显示 CMORPH 卫星 数据 对 比 实 
况 观测 ,其 降水 范围 显示 出 明显 的 高 佑 现象 ,尤其 
在 新 疆 北部 的 阿勒泰 . 塔 城 以 及 新 疆 南 部 的 喀什 、 
和 田地 区 ,其 融合 产品 CMPA 在 新 疆 北部 地 区 有 了 
很 大 改进 , 而 在 新 疆 南部 地 区 仍 存在 大 量 高 估 区 
域 ,而 图 8c 融合 降水 pdf-oi_CRT 则 在 2 个 区 域 都 做 
出 了 很 好 地 处理 ,与 实况 降水 空间 分 布 较为 一 致 。 
由 伊犁 河谷 降水 中 心 来 看 ,CMPA 也 比 实况 明显 偏 
大 ,由 于 CMORPH 卫星 降水 在 新 疆 地 区 存在 明显 的 
高 估 ™, 因 此 ,本 研究 基于 IMERG 的 融合 试验 有 望 


得 到 更 为 准确 的 区 域 降水 量 。 
4 结论 


(1) 概率 密度 匹配 法 (PDF) 对 原始 IMERG 的 系 
统 误差 有 很 好 的 订正 效果 ,使 原始 卫星 降水 产品 的 
绝对 平均 误差 与 相对 误差 均 明 显 减 小 。 

(2) 交叉 验证 结果 显示 ,在 IMERG 降水 产品 系 
统 误差 经 PDF 订 正 后 显著 降低 的 基础 上 ,进一步 利 
用 0I 方 法 与 实况 降水 融合 ,使 最 终 得 到 的 格 点 降水 
数据 集 质量 进一步 提高 。 

(3) 研究 体现 了 PDF-0I 方 法 应 用 于 IMERG T 
星 降水 产品 的 有 效 性 , 且 融 合 产 品 的 降水 量 与 CM- 
PA 融合 降水 相 比 与 实况 更 为 一 致 。 
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Multisource precipitation data merging experiment in Xinjiang 


LU Xin-yu, LIU Yan', WANG Xiu-qin', LIU Jing, WEI Ming’, 
SONG Zhi-guo, ZHANG Ying-xin' 

(1. Institute of Desert Meteorology, China Meteorological Administration, Urumqi 830002, Xinjiang, China; 
2. Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing 
University of Information Science and Technology, Nanjing 210044, Jiangsu, China; 

3. Yishui meteorological bureau of Shandong, Linyi 276400, Shandong, China) 


Abstract: Precipitation is a fundamental component of the Earth’ s water cycle. Recognizing the temporal and 
spatial distributions of precipitation in a specific area is helpful to better understand the ecological and 
environmental changes in this area, and the most important of these is to accurately obtain the precipitation data for 
the entire area. Therefore, developing high spatiotemporal resolution and high-quality precipitation data is an urgent 
need as the development of modern meteorological services. In this paper, the study on a fusion method based on 
two-source precipitation, that is, satellite precipitation data and ground observation, is carried out in the Xinjiang 
area with sparse stations in China, which can be of great significance to accurately reflect the precipitation 
distribution information and the characteristics of precipitation evolution process in this area. 

Using hourly precipitation data in summer (June- August) 2016 of about 2000 automatic weather stations (AWS) 
over Xinjiang and using the latest generation of GPM IMERG satellite precipitation data as the first guess, a two- 
step strategy based on the probability density function (PDF) matching and optimal interpolation (OI) technique 
(PDF-OD was conducted to a study on fusion methods of precipitation in arid regions. First, the bias of the IMERG 
precipitation is corrected by matching the PDF of the observation with the satellite data. Then, the bias-corrected 
IMERG precipitation (pdf_CRT) is used as the initial estimation field and the observation field (OBS) as the 
truthvalue; the combination of the pdf_CRT and OBS through the OI method is carried out to finally obtain hourly 
areal precipitation data (pdf-oi_CRT) in Xinjiang. The cross-validation results show that, compared with the OBS 
and the satellite precipitation, the PDF-OI fusion method not only greatly eliminates the bias but also significantly 
improves the accuracy of the original satellite precipitation data. 

In this study, a useful exploration is made in developing a grid precipitation data set in Xinjiang, which is 
located in a middle-to high-latitude arid area. The results show that the study on the multisource precipitation fusion 
method helps develop grid precipitation data set. The high temporal and spatial resolution grid precipitation data set 
will help us understand the characteristics and occurrence and development regularity of hourly heavy precipitation 
in the Xinjiang area. Furthermore, it will also provide data support for improving the accuracy of forecast and 
preventing secondary meteorological disasters and the technical foundation for developing high spatial and temporal 
resolution precipitation products in the Central Asian region. 


Keywords: satellite-gauge merging; hourly precipitation; fusion; Xinjiang 


